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Abstract 

Using a sample of 3 x 10 6 ip(2S) decays recorded by the CLEO detector, we study three body 
decays of the XcO, Xci> an d Xc2 produced in radiative decays of the 4>{2S). We consider the decay 
modes ir + iT~~ri, K + K~n, ppn, 7r + 7r~?/, K + K~it°, ppir , ir+K^K®, and K + pk measuring branching 
fractions or placing upper limits. For Xd — * ir + ir~r], K + K~ir°, and tt + K~ K§ our observed samples 
are large enough to study the substructure in a Dalitz plot analysis. The results presented in this 
document are preliminary. 
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Decays of the Xcj states are not as well studied both experimentally and theoretically as 
those of other charmonium states. It is possible that the color-octet mechanism, ccg — > 2(qq) 
could have large effects on the observed decay pattern of the Xcj states Thus any 
knowledge of any hadronic decay channels for these state is valuable. 

CLEO has gathered a large sample of e + e~ — > i>(2S) which leads to copious production 
of the Xcj states in radiative decays of the if)(2S). This contribution describes our study of 
selected three body hadronic decay modes of the Xcj to two charged and one neutral hadron. 
This is not an exhaustive study of Xcj hadronic decays; we do not even comprehensively 
cover all possible h + h~h° decays, but simply take a first look at the rich structure of Xcj 
decays in our initial ip(2S) data sets. With the CLEO III detector configuration 2j], we have 
observed an integrated luminosity of 2.57 pb _1 and the number of '0(25') events is 1.56 x 10 6 . 
With the CLEO-c detector configuration 3] we have observed 2.89 pb _1 , and the number 
of events is 1.52 x 10 6 . Note that the apparent mis-match of luminosities and event totals 
is due to different beam energy spreads for the two data sets. 

Our basic technique is an exclusive whole event analysis. A photon candidate is combined 
with three hadrons and the 4-momentum sum constrained to the known beam energy and 
small beam momentum caused by the beam crossing angle taking into account the measured 
errors on the reconstructed charged tracks, neutral hadron, and transition photon. We 
cut on the x 2 °f this fit, which has four degrees of freedom, as it strongly discriminates 
between background and signal. Efficiencies and backgrounds are studied in a GEANT- 
based simulation p[ of the detector response to underlying e + e~ — > ip{2S) events. Our 
simulated sample is roughly ten times our data sample. The simulation is generated with a 
1 + A cos 2 9 distribution in cos 9, where 9 is the radiated photon angle relative to the positron 
beam axis. A El transition, as expected for i/j(2S) — > JXcj, implies A = 1, —1/3, +1/13 
for J = 0,1,2 particles. The efficiencies we quote use this simulation. The differences of 
efficiencies due to various 9 distributions are negligible as we accept transition photons down 
to our detection limit. 

A kinematic fit is made to each event. For most modes we select events with an event 4- 
momentum fit x 2 l ess than 25, but background from ip(2S) — > J/^7r°7r° followed by charged 
two body decays of the J/ip with one of the 7r° decay photons lost fakes ip(2S) — > •"fXcj 
'jppTT . For this mode the cut on x 2 is tightened to 12. The efficiency of this cut is ~ 95% 
for all modes except ppn° where it is ~ 80%. Little background survives. 

Photon candidates are energy depositions in our crystal calorimeter that have a transverse 
shape consistent with that expected for an electromagnetic shower without a charged track 
pointing toward it. They have an energy of at least 30 MeV. Transition photon candidates 
are vetoed if they form a tt° or rj candidate when paired with a second photon candidate. 
Photon candidates that are formed into neutral particles further must have an energy of 
more than 50 MeV if they are not in the central barrel of our calorimeter, it — > 77 and 
rj — > 77 candidates are formed from two photon candidates that are kinematically fit to the 
known resonance masses using the event vertex position, determined using charged tracks 
constrained to the beam spot. We select those giving a x 2 from the kinematic mass fit with 
one degree of freedom of less than 10. 

We also select the rj — > 7r + 7r~7r mode combining two charged pions with a 7r° candidate 
increasing the efficiency for 7] reconstruction by about 25%. The same sort of kinematic mass 
fit as used for 7r°'s and rj — > 77 is applied to this mode, and again we select those giving 
a x 2 °f l ess than 10. Similarly we combine the mass-constrained 77's together with two 
charged pions to make rj' candidates, mass constrain them, and select those with x 2 < 10. 
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In addition, we include the decay mode rf — > 7p. Here the background is potentially high 
because of the large number of noise photons, so we require ^photon > 200 MeV. In addition, 
we require the 7r + n~ mass to be within 100 MeV/c 2 of the p mass. 

Charged tracks have standard requirements that they be of good fit quality. Those coming 
from the origin must have an impact parameter with respect to the beam spot less than the 
maximum of (5.0 — (3.8 • p)) mm or 1.2 mm, where p is the measured track momentum in 
GeV/c. Kg — > 7r + 7r~ and A — > pn~ candidates are formed from good quality tracks that are 
constrained to come from a common vertex. The K s flight path is required to be greater 
than 5 mm and the A flight path greater than 3 mm. The mass cut around the K s mass is 
±10 MeV/c 2 , and around the A mass ±5 MeV/c 2 . Events with only the exact number of 
selected tracks are accepted. This selection is very efficient, >99.9%, for events passing all 
other requirements. 

Pions are required to have specific ionization, dE/dx, in our main drift chamber within 
four standard deviations of the expected value for a real pion at the measured momentum. 
For kaons and protons, a combined dE/dx and RICH likelihood is formed and kaons are 
simply required to be more kaon-like than pion or proton-like, and similarly for protons. 
Cross feed between hadron species is negligible after all other requirements. 

In modes comprising only two charged particles, there are some extra cuts to eliminate 
QED background which produce charged leptons in the final state. Events are rejected if 
the sum over all the charged tracks produces a penetration into the muon system of more 
than five interaction lengths. Events are rejected if any track has 0.92 < E/p < 1.05 and it 
has a dE/dx consistent with an electron. This latter cut is not used for pp modes because 
anti-protons tend to deposit all their energy in the calorimeter. These cuts are essentially 
100% efficient for the signal, and do eliminate the small QED background. 

The efficiencies averaged over the CLEO III and CLEO-c data sets for each mode includ- 
ing the branching fractions rj — > 77, 77 — > 7r + 7r _ 7r°, and rf — > i]7t + tt~ are given in Tables ITl lIIII 
for Xco, Xci, and Xc2 respectively. Figures [T]|H] show the mass distributions for the eight 
Xcj decay modes selected by the analysis described above. Signals are evident in all three 
Xcj states, but not in all the modes. Backgrounds are small. The mass distributions are 
fit to three signal shapes, Breit-Wigners convolved with Gaussian detector resolutions, and 
a linear background. The Xcj masses and intrinsic widths are fixed at the values from the 
Particle Data Group compilation jj^]. The detector resolution is taken from the simulation 
discussed above. The simulation properly takes into account the amount of data in the two 
detector configurations, and the distribution of different decay modes we have observed. We 
approximate the resolution with a single Gaussian distribution, and variations are consid- 
ered in the determination of systematic uncertainty. The detector resolution dominates for 
the Xci and Xc2, but is similar to the intrinsic width of the Xco- The fits are displayed in 
Figures and summarized in Tables Hi ll I II Note that for the Xco i n Table U the five modes 
for which no significant signal are found are forbidden by parity conservation. 

We consider various sources of systematic uncertainties on the yields. We varied the 
fitting procedure by allowing the Xcj masses and intrinsic widths to float. The fitted masses 
and widths agree with the fixed values from the particle data group, and we take the maxi- 
mum variation in the observed yields, 4%, as a systematic uncertainty from the fit procedure. 
Allowing a curvature term to the background has a negligible effect. For modes with large 
yields we can break up the sample into CLEO III and CLEO-c data sets, and fit with reso- 
lutions and efficiencies appropriate for the individual data sets. We note that the separate 
data sets give consistent efficiency corrected yields and the summed yield differs by 2% from 
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FIG. 1: Mass distribution for candidate \cj — * tt + tt r]. The displayed fit is described in the text. 



TABLE I: Parameters used and results of fits to the Xcj mass distributions of Figures UM\ for the 
XcO signal. The fit is described in the text and the yield error is only statistical. If no significant 
signal is observed the 90% confidence level limit is also shown. 



Mode 


Efficiency (%) 


Resolution (MeV) 


Yield 


K+K-j] 
pprj 

■K + ir~r]' 
K+R-ir 

ppTT 

K+pA 


17.6 
13.6 
15.8 
7.8 
26.4 
27.1 
19.0 
16.7 


6.23 
6.10 
5.42 
4.38 
6.47 
5.80 
4.75 
4.38 


4.0 ±4.1 (< 10.6) 

3.1 ± 2.7 (< 9.3) 

i 7 7+5.2 
i ' -'-4.8 

1.8 ±3.8 (< 8.5) 
-3.3 ±2.7 (< 4.7) 
46.2±?;g 

0.0 ± 1.0 (< 2.7) 
51.2t 8 7 i 
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FIG. 2: Mass distribution for candidate Xcj — ► K + K rj. The displayed fit is described in the text. 



TABLE II: Parameters used and results of fits to the XcJ mass distributions of Figures WffiQ for the 
Xci signal. The fit is described in the text and the yield error is only statistical. If no significant 
signal is observed the 90% confidence level limit is also shown. 



Mode 


Efficiency (%) 


Resolution (MeV) 


Yield 




18.1 


5.85 


255+Iq 


K+K-r] 


14.2 


5.66 


i4.i±!:i 


pprj 


17.6 


5.41 


3.2 ±2.3 (< 7.6) 




8.3 


4.37 


57.61$ 


K+K-tt 


28.2 


6.79 


157 ±13 


PpTT° 


29.5 


5.23 


q q+ 3 - 8 
y - y -3.2 


tt+K~K^ 


19.8 


4.37 


249 ± 16 


K + pA 


17.5 


4.38 


16 3+ 47 
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FIG. 3: Mass distribution for candidate Xcj —> PPV- The displayed fit is described in the text. 



TABLE III: Parameters used and results of fits to the \cj mass distributions of Figures for the 
Xc2 signal. The fit is described in the text and the yield error is only statistical. If no significant 
signal is observed the 90% confidence level limit is also shown. 



Mode 


Efficiency (%) 


Resolution (MeV) 


Yield 




17.8 


5.58 


26.211$ 


K+K-j] 


14.1 


5.53 


6.9 ±2.9 (< 12.54) 


pprj 


16.9 


5.08 


q 5+3.8 
y -°-3.0 


7r + 7r~ 7/ 


8.2 


4.32 


12.4±tf 


K+K-7T 


27.5 


6.85 


24.8±|-f 




28.9 


5.10 


37.1+1? 




17.2 


4.45 


36.8±g;g 


K+pA 


17.5 


4.32 
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FIG. 4: Mass distribution for candidate Xcj - > tt + tt r( . The displayed fit is described in the text 



the standard procedure, which is small compared to the 8% statistical uncertainty. We take 
this as the systematic uncertainty from our resolution model. From studies of other pro- 
cesses we assign a 0.7% uncertainty for the efficiency of finding each charged track, 2.0% for 
the 77 resonances, 1.0% for each extra photon, 1.3% for the particle identification for each 
K and p, 2.0% for secondary vertex finding, and 3.0% from the statistical uncertainty on the 
efficiency determined from the simulation. We study the cut on the x 2 of event 4-momentum 
kinematic fit in the three large yield Xci signals by removing the \ 2 cut, selecting events 
around the xa mass peak, subtracting a low mass side band, the only one available, and 
comparing the simulated \ 2 distribution for signal events with the data distribution. This 
comparison is shown in Figured The agreement between the data and simulation is good, 
and comparing the inefficiency introduced by our cut on the 4-momentum kinematic fit x 2 
between the data and the simulation we assign a 3.5% uncertainty on the efficiency due 
to uncertainty in modeling this x 2 distribution. The simulation was generated assuming 
3-body phase-space for the Xcj decay products. Deviations from this are to be expected. 
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K + K 7T°. The displayed fit is described in the 



Based on the results of the Dalitz plot analyses discussed below we correct the efficiency 
in the xa — * 7r + 7r~?7, Xci —* K + K~tt°, and xa — > K~7i + Kg modes by a relative +6% 
to account for the change in the efficiency caused by the deviation from a uniform phase 
space distribution of decay products to what we actually observe. This correction only has 
a noticeable impact on the Xci K^ir + Kg mode. 

To calculate Xcj branching fractions, we use previous CLEO measurements for the 
ip(2S) -> Xcj branching fractions of 9.3±0.14±0.61%, 9.07±.11±.54% and 9.22±.11±.54% 
for J— 0,1,2 respectively 0|. The uncertainties on these branching fractions are included in 
the systematic uncertainty on the Xcj branching fractions we report. Preliminary results for 
the three body branching fractions are shown in Table IIV1 Where the yields do not show 
clear signals we calculate 90% confidence level upper limits using the yield central values 
with the statistical errors from the yield fits combined in quadrature with the systematic 
uncertainties on the efficiencies and other branching fractions. We assume the uncertainty 
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FIG. 6: Mass distribution for candidate Xcj —> ppn ■ The displayed fit is described in the text. 



TABLE IV: Preliminary branching fractions in %. Uncertainties are statistical, systematic due to 
detector effects plus analysis methods and a separate systematic due to uncertainties in the ip(2S) 
branching fractions. Limits are at the 90% confidence level. 



Mode 


XcO 


Xcl 


Xc2 


ir + ir~r] 




< 0.021 


0.52 ± .03 ± .03 ± .03 


0.051 ±.011 ±.004 ±.003 


K+K-T] 




< 0.024 


0.034 ± .010 ± .003 ± .002 


< 0.033 


PpT] 


0.038 ± 


.010 ± .003 ± .02 


< 0.015 


.019 ± .007 ± .002 ± .002 






< 0.038 


0.24 ± .03 ± .02 ± .02 


< 0.053 






< 0.006 


0.200 ± .015 ± .018 ± .014 


0.032 ± .007 ± .002 ± .002 


ppTT 


0.059 ± 


010 ± .006 ± .004 


0.014 ± 0.005 ± 0.001 ± 0.001 


0.045 ± .007 ± 0.004 ± .003 


ir+K-K 




< 0.010 


0.84 ± .05 ± .06 ± .05 


0.15 ±.02 ±.01 ±.01 


K + pA 


0.114 ± 


016 ± .009 ± .007 


0.034 ± .009 ± .003 ± .002 


0.088 ± .014 ± .07 ± .006 
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FIG. 7: Mass distribution for candidate Xcj KgK tt + . The displayed fit is described in the 
text. 



is distributed as a Gaussian and the upper limit is the branching fraction value at which 
90% of the integrated area of the Gaussian falls below. We exclude the unphysical region, 
negative branching fractions, for this upper limit calculation. We note that the ratio of rates 
expected from isospin symmetry, as discussed in Appendix HI Equations fTol and 123} expected 
to be 4.0 is consistent with our measurement: 

rtxa^ifiK+K-) -4-2 ±0.7. (1) 

We choose the three high-statistics signals xa —> ^n^r], Xci —* K + K~tt°, and x& ~^ 
K~-K + Kg for Dalitz plot analysis to study resonance substructure. For the Dalitz analysis 
only those events within 10 MeV, roughly two standard deviations, of the signal peak are 
accepted. For Xci vr + 7r _ ?7 there are 228 events in this region and the signal fit finds 224.2 
signal events and 5.1 combinatorial background. For Xci K + K~tt° there are 137 events 
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FIG. 8: Mass distribution for candidate Xcj K + pA. The displayed fit is described in the text. 
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FIG. 9: Distribution of the x 2 °f the beam energy constrained mass fit. This is shown with a mass 
cut around the Xd signal region and after a data sideband subtraction. The plot on the left is for 
the 7T + 7r~77 mode, center is the K + K~tt° mode, and right is the 7t + K~Kq i mode. The data are 
shown by points and the simulation of signal events is shown by the histogram. 
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accepted with the fit finding 137.8 signal and 2.4 background events, and for Xci — > K~tt + K S} 
the numbers are 234 events, of which 233.2 are signal and 0.8 are background. In all cases 
the contribution from the tail of the Xc2 is less than one event. 

An unbinned maximum likelihood fit is used along with other methods as described 
in 0] in order to perform the Dalitz plot analysis. We only summarize our methods here. 
Efficiencies are determined with simulated event samples generated uniformly in phase space, 
and run through the analysis procedure described above. The efficiency across the Dalitz 
plots is fit to a two dimensional polynomial of third order in the Dalitz plot variables. The 
fits are of good quality and the efficiency is generally flat across the Dalitz plot. 

When fitting the data Dalitz plot the small contributions from backgrounds are neglected. 
We use an isobar model to describe resonance contributions to the Dalitz plots taking into 
account spin and width dependent effects. Narrow resonances are described with a Breit- 
Wigner amplitude with the resonance parameters taken from previous experiments 0. For 
the scalar resonances a (980) and / (980) we use a Flatte parameterization. We use the 
ao (980) line-shape from the Crystal Barrel Collaboration [8| and the details of the /o(980) 
are unimportant as it is used only in systematic studies. For low mass 7r + 7r~(o~) and K7t(k) 
S-wave contributions we choose a simple description, one which is adequate for our small 
sample [^]. 

We are examining the e + e~ — > ip{2S) — > JXa process. In such a decay the xa should be 
polarized. In principal a more complete analysis would take into account the angle of the 
photon with respect to the e + e~ beams collision axis and decompose the Xci decay into its 
partial waves. We barely have the statistics to do a reasonable Dalitz analysis suggesting 
that a higher dimensional partial wave analysis would be hopeless. See Appendix UJJ for a 
discussion and the formalism of how polarization would affect the Dalitz analysis when the 
intermediate resonance is not spin zero even when integrated over the random polarization 
direction of the xa- Keeping this in mind, for this Dalitz plot analysis we decided to use 
angular distributions from 10]. We have tested different angular distributions and include 
the variations clS db systematic uncertainty. 

Figure ITU1 shows the Dalitz plot and three projections for Xci —> 7r + 7r - ?7. There are clear 
contributions from a (980) ± 7r T and / 2 (1270)?7 intermediate states, and significant accumu- 
lation at low 7r + 7r~ mass. Note that the a (980) can contribute in two decay modes to the 
Dalitz plot. An isospin Clebsch-Gordan decomposition for this decay, given in Appendix UJ 
Equation EE shows that amplitudes and strong phases of both charge-conjugated states 
should be equal. The overall amplitude normalization and one phase are arbitrary parame- 
ters and we set a ao{980) + = a ao{980) - = 1, ao ( 98 o)+ = 0a o (98O)- = 0. All other fit components 
are defined with respect to these choices for ao(980). 

Our initial fit to this mode includes only ao(980) ± 7r =F and f2(l270)i] contributions, but 
has a low probability of describing the data, 0.13%, due to the accumulation of events at 
low 7r + 7r~ mass. To account for this we try Kg, / (980), p(770), and a resonances. Only the 
p(770) and the a give high fit probability, 49% and 58% respectively. However the decay 
Xci p(770)?7 is C-forbidden, and the low mass distribution is not well represented by p(770) 
which only gives an acceptable fit due to its large width and the limited statistics of our 
sample. The a describes well the low tt + tt~ mass spectrum, and we describe the Dalitz plot 
with ao(980) ± 7T =F , /2(1270)?y, and at] contributions. Table IVl gives the preliminary results 
of this fit which has a probability to match the data of 58.1%. Variations to this nominal 
fit give the systematic uncertainties shown in the table. We allow the 2D-efficiency to vary 
with its polynomial coefficients constrained by the results of the fit to the simulated events; 
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FIG. 10: Dalitz plot and projections on the three mass squared combinations for Xd vr + 7r r\. 
The displayed fit projections are described in the text. 



TABLE V: Preliminary fit results for Xd i]tt + tt Dalitz plot analysis. The uncertainties are 
statistical and systematic. 



Contribution 


Amplitude 


Phase (°) 


Fit Fraction (%) 


a (980) ± 7r =F 


1 





56.2 ±3.6 ±1.4 


/ 2 (1270)i7 


0.186 ±0.017 ±0.003 


-118 ± 10 ±4 


35.1 ±2.9 ± 1.8 


arj 


0.68 ±0.07 ±0.05 


-85 ±18 ±15 


21.7 ±3.3 ±0.5 
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the mass of the ao(980) and its coupling constants are allowed to float, the parameters of the 
cr-pole are allowed to float, and we allow additional contributions from p(770)?7, /o(980)r/, 
KgT], and 7Ti(1400)?7 . The deviation from the nominal fit over these variants gives the 
systematic uncertainties shown in Table For the additional contributions we do not 
observe amplitudes that are significant and we limit their fit fractions at the 95% confidence 
level to p(770)r] < 12%, / (980)7/ < 3.5%, K° s r) < 0.8%, 7Ti(1400)t7 < 3.2%. We note that 
with higher statistics this mode may offer one of the best measurements of the parameters 
of the a (980). 

The Dalitz plot for \ci — * K + K~tt° decay and its projections are shown in Figure fTTl 
and for xa —> n + K~Kg in Figure We do a combined Dalitz plot analysis to these 
modes taking advantage of isospin symmetry. An isospin Clebsch-Gordan decomposition 
for these decays, Appendix HI Equations [TTJ [T21 and EH shows that these two Dalitz plots 
should have the same set of amplitude and phase parameters for all K*K and a (980)7r 
intermediate states. The relative factor — v^2 between two Dalitz plot amplitudes does 
not matter due to the individual normalization of their probability density functions. In 
the combined fit to these two Dalitz plots, we use the following constraints on amplitudes 
and phases: a^*+ = a^*- = clk*° — oq^o = 4>k*+ — 4>k*- — 4>k*° — 4>j^o = 4>k*, 

Oa(980)+ = a a (980)- = O-a(980)° = «-a(980) , and a (98O)+ = 0a(98O)~ = 0a(98O)° = 0a(98O)- The Over- 
all amplitude normalization and one phase are arbitrary parameters and we set ax* (892) — 1 
and 0k*(892) = 0. 

The limited size of this sample, even in the combined Dalitz plot analysis, and the 
many possible contributing resonances leave us unable to draw clear conclusions. Visual 
inspection shows apparent contributions from K* (892) ± K T , K*(892)°K° S , J R'*(1430) ± ir =F , 
K*(U30)°K° S , a (980) 7r°, and a (980) ± vr T . It is not clear if the K*(U30) are K* or K*, 
and many other Kir and KK resonances can possibly contribute. Our best fit preliminary 
result is shown in Table EU showing statistical and systematic errors. This fit has a good 

TABLE VI: Preliminary results of the combined fits to the Xd ~^ K + K~tt° and Xci KgKir 
Dalitz plots. 



Contribution 


Amplitude 


Phase (°) 


Fit Fraction (%) 


K*(892)K 
K £(1430) K 

k{(uw)k 

K*(lQm)K 
a (980)vr 


1 

0.50 ±0.09 ±0.12 
5.3 ± 1.0 ±0.1 
2.3 ±0.5 ±0.5 
10.8 ± 1.2 ± 1.2 




-2 ±13 ±6 
77 ± 12 ± 16 
-38 ± 12 ± 12 
-112 ± 12 ±3 


19.7 ±4.0 ±2.0 
18.0 ±6.6 ±3.2 
36.0 ± 12.8 ±3.0 
11.2 ±5.4 ±2.7 
29.5 ±7.1 ±2.6 



probability of matching the data, 32.1%, and agrees with fits done to the separate Dalitz 
plots not taking advantage of isospin symmetry. Unfortunately we can change our decay 
model by swapping the i^*(1430) plus if* (1680) contributions for a flat non-resonant back- 
ground, or a k and get fits of lower, but still acceptable, probability of matching the data. 
We also find that an alternative solution using the same set of contributions fits the data 
acceptably, but less well than the displayed result. This alternative result disagrees with the 
nominal result by more than the statistical uncertainties should allow. We conclude that 
our Xci — * KKtt sample is too small to extract the resonant substructure. We do clearly see 
contributions from K*(892)K and ao(980)7r at roughly the 20% and 30% level respectively. 
The balance of the Dalitz plots are probably from higher mass resonances. 

In conclusion we have searched for and studied selected three body hadronic decays of the 
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FIG. 11: Dalitz plot and projections on the three mass squared combinations for Xci K + K 7r c 
The displayed fit projections are described in the text. 



Xco, Xci, Xc2 produced in radiative decays of the ip(2S) in e + e~ collisions observed with the 
CLEO detector. Our preliminary observations and branching fraction limits are summarized 
in Table HV1 In Xci vr + 7r _ r] we have studied the resonant substructure using a Dalitz plot 
analysis, and our preliminary results are summarized in Table Similarly in Xci —* KKir 
we clearly see contributions from K*(892)K and ao(980)7r at roughly the 20% and 30% level 
respectively. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. D. Cronin-Hennessy and A. Ryd thank the A. P. Sloan 
Foundation. This work was supported by the National Science Foundation, the U.S. De- 
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The displayed fit projections are described in the text. 
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I. APPENDIX: CLEBSCH-GORDAN DECOMPOSITION FOR X 



cl 



In order to constrain amplitudes and phases in Xci decays we use a Clebsch-Gordan 
decomposition of the Xci (the state with \I = 0, Iz = 0)) for possible isospin subsystems: 

Xd - (Kn) I=l/2 K, (2) 

Xd - (Ktt) i=1/2 K, (3) 
Xci - {KK) I=1 ir. (4) 
Below we use Clebsch-Gordan decomposition rules, \ J,M) = J2f c f\ m i, %)/ from Ref. 



A. Clebsch-Gordan decomposition for K* — » Kir decays 



We assume that K*s with 1=1/2 form two isodoublets: (K* + , K*°) and (K*°, K*~) with 
{Iz = \ jz = — |) respectively. The Clebsch-Gordan decomposition rules for isospin states 
of product particles 1 X 2 £1X6! 



II I) = Jl\\ -1) _ A /i |o i) 

11 _I\ = ./I|o —I) _ -1 1) 

I 2' 2/ V 3 1 ' 2/ V 3 1 '2/ 
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B. Cases of Xd -» (Kir) I=1 / 2 K and %ci -> (Kn) I=1 / 2 K decays 

For Xd — > and Xci — > if*if modes we use | x | rule: |0, 0) = ^(||, — f) 



.5» 



l l 

2' 2 



Xcl 



Xcl 



1 

1 

71 



(K* + K- - K*°K°), 
(K*K° - K*-K + ). 



Combining Equations 191101 with Equations OBI we get 



Xcl 



V2 



X, 



■iV2 



Assuming charge symmetry the amplitudes in Equations ITT1 and H"2l should be equal, 
these equations we may get the ratio of rates: 

r( Xc i - n + K~K°)/T( X ci - n°K + K-) = 2, 

r( X ci -> 7r-^+^)/r( Xcl -> 7r°K+i<:-) = 2, 



or their sum 



r(x 



cl — > 7T 



~K~K°) + r(x c i -> n-K + K°) = 4T( Xc i -> A+iT 



(9) 
(10) 

(11) 

(12) 
From 

(13) 
(14) 

(15) 



C. Case of Xcl — ► {KK)j = itt decay 

For Xci - ► 07r modes we use the lxl rule: |0, 0) = — 1) — |0, 0) + | — 1, 1)) 



Xcl 



' + - i - - 

a 7r — a n + a ir 



For a — > If if we use the 2 x 2 rmes: 
1 i\ 



a+: |1,1) = || i) 



a°: |l,0) = ^(|l,-i) + 



a-. |1,-1) = 



- 1 !)) 

2' 2// 



v/2 



(K + K~ + K°K°), 



Combining Equation ITol with Equations IT71ITU1 we get 

Xciv^ -> (if+Zo)7r- - -L [(K + if-)vr° + (K°K°)n° 

V2 . 



(16) 

(17) 

(18) 
(19) 

(20) 
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2, (21) 
2, (22) 

-> n°K + K-). (23) 
D. Consequences for Dalitz plot analysis 

Comparing Equations IT3llTo1 for intermediate states with K* and Equations l2~Tll2*3l for 

intermediate states with a(980), one may notice that they are identical. Thus observations 
of 7T + K~K° and charge conjugated tt~K + K° final states on the same Dalitz plot will yield 
the same ratio between all K* and a(980) amplitudes for both K°K~n + and K + K~tt° 
Dalitz plots. The ratio of amplitudes between these two Dalitz plots is —y/2. The relative 
negative sign between amplitudes does not matter, because the rates depend on the matrix 
element squared. The relative factor between amplitudes does not matter, because the 
normalizations are different for each Dalitz plot. This isospin analysis implies that these two 
Dalitz plots, K°K~tt + and K + K~tt°, can be parametrized using a common set of parameters 
for each of K*s and a(980) intermediate state with 

a K *+ = a K *- = a K *o = a-^za = a K * } (24) 

4>K*+ = 4>K*~ = 4>K*o = = <j>K*, (25) 

a a(980)+ — a a(980)" = a a(980)° = «a(980); (26) 
0a(98O)+ = 0a(98O)" = 0a(98O)° = 0a(98O)- (27) 



From Equations |20] we may get the ratio of rates: 

r( Xc i -> 7r+K-K°)/r(xci - n°K+K-) = 
r(xci -> 7i-K + W)/T(xci -> n°K + K-) = 

or their sum 

r(xci ^ n + K-K°) + T( Xc i ^ n-K+W) = 4T( X d 



II. APPENDIX: THE EFFECT OF X d POLARIZATION 
A. Angular distributions 

In this analysis we use the angular distributions explicitly shown in Table IVTD derived for 
a non-polarized decaying particle (x&) m R- e f- 0|- The parameters in Table IVTT1 follow the 
conventions of the original publication: J - spin of the initial particle J — > j + c; j - spin 
of the resonance j — > a + b, where a,b,c - spin particles; L - orbital angular momentum 
between (j p ) resonance and the recoil particle c. A relativistic correction factor z and cos 8, 
where 9 is an angle between directions of particles a and c in resonance j p rest frame is 
discussed in Appendix III Bl 

In order to account for the polarization of the initial Xci produced in the radiative decay 
ip(2S) — > 7X c i, a full Partial Wave Analysis formalism, for example from Ref. 0], would 
be needed. The matrix element amplitudes would depend on additional angular variables 
in addition to the two invariant squared masses as in the regular Dalitz plot analysis. The 
number of parameters would also be increased to account for different helicity amplitudes 
contributions. On the other hand, an expected X & production angular distribution, oc 
1 — | cos 2 9 Xcl , is ~15% different from uniform and we expect that the effect of polarization 
is small in the decays under study. 
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TABLE VII: Angular distributions. 



J - 


> 7 +£ 


Probability 


0- 


+0+0 


uniform 


1- 


+0+1 


uniform 


2- 


+0+2 


uniform 


0- 


-+1+1 


(l + z 2 )cos 2 6> 


0- 


+2+2 


(z 2 + 3/2) 2 (cos 2 6»- 1/3) 2 


1- 


+1+0 


1 + z 2 cos 2 9 


1- 


+1+1 


sin 2 6 


1- 


+1+2 


1 + (3 + 4z 2 ) cos 2 


1- 


+2+1 


(1 + z 2 )[l + 3 cos 2 9 + 9z 2 (cos 2 9 - 1/3) 2 ] 


1- 


+2+2 


(l + .z 2 )cos 2 0sm 2 (*) 


2- 


+1+1 


3 + (l + 4z 2 )cos 2 6l 


2- 


+1+2 


sin 2 9 


2- 


+2+0 


1 + z 2 /3 + z 2 cos 2 9 + z 4 (cos 2 9 - 1/3) 2 


2- 


+2+1 


1 + z 2 /9 + (z 2 /3 - 1) cos 2 9 - z 2 (cos 2 9 - 1/3) 2 


2- 


+2+2 


1 + z 2 /9 + (z 2 /3 - 1) cos 2 9 + (16z 4 + 21z 2 + 9)(cos 2 9 - l/3) 2 /3 (*) 



(*) These formulas have been derived based on the covariant helicity formalism approach also 

discussed in Ref. [lo| . 

B. Expressions for relativistically non-invariant variables in terms of Dalitz plot 
invariant variables 



Although a covariant spin tensor formalism is applied, for simplicity the formulas in 
Ref. H3 and Table |VlTl are expressed in terms of relativistically non-invariant variables z 
and cos 9. Here we discuss the meaning of these variables and their expression in terms of 
particle masses and invariant masses. 

Assuming the decay d — > Rc, R -+ ab we may derive kinematic variables used in for- 
mulas for angular distribution in terms of Dalitz plot invariant variables. The energy and 
momentum of particle a in the resonance R or (ab) rest frame is a trivial expression 



E n 



™ab + m a 



2m n 



Pn 



E 2 



(28) 



The energy and momentum of particle c in the resonance R rest frame can be obtained 
from the invariant mass m 2 d expression in the R rest frame. Indeed, m 2 = (p c + pn) 2 , in the 



resonance R rest frame p c = (E C} P C ) } p R 



E r 



m 



ab 



m 



(m R ,0), 

2 



m 2 + mil + 2mnE c , hence 



2m, 



Pr 



ab 



E 2 — m 2 . 



(29) 



Now, the cos 9 between directions of particles a and c in resonance R rest frame can be 
expressed through the known energies and momenta of particles a and c and their measured 
invariant mass squared m 2 



l ac- 



cos 9 



m 2 a + m 2 c + 2E a E c 

2P n P r 



(30) 
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The resonance energy and relativistic correction factor (used in in Table fVTT]) in the decaying 
particle d rest frame are: 

m g +m ; t - m ; , 2 = 7| _ 1 = ^_ L (31) 

2m d n m 2 ab 
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